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Abstract 

We extend the study of a model of competitive cluster growth in an active medium from a 
regular topology to a complex network topology; this is done by adding nonlocal connections with 
probability p to sites on a regular lattice, thus enabling one to interpolate between regularity 
and full randomness. The model on networks demonstrates high sensitivity to small changes in 
initial configurations, which we characterize using damage spreading. The main focus of this 
paper is, however, the devising of survival strategies through selective networking, to alter the fate 
of an arbitrarily chosen cluster: whether this be to revive a dying cluster to life, or to make a 
weak survivor into a stronger one. Although such goals are typically achieved by networking with 
relatively small clusters, our results suggest that it ought to be possible also to network successfully 
with peers and larger clusters. The main indication of this comes from the probability distributions 
of mass differences between survivors and their immediate neighbours, which show an interesting 
universality; they suggest strategies for winning against the odds. 
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I. INTRODUCTION 



The study of complex systems with interacting agents is of great interest in very many 
fields ranging from physics through quantitative biology to economics [Ij. Modelling the 
dynamics of interacting species should capture the competition between the species for 
growth and survival, which underlies all evolutionary dynamics. One of the most famous 
models of this class is the Lotka-Volterra model used for studying predator-prey interactions 
[2!, which is loosely predicated on the principle of the survival of the fittest [s]. 

Typically, little attention has been paid to the case where agents interact through an 
active medium; this is a case where the Lotka-Volterra model cannot be simply applied. A 
very recent study [4] on microbial populations strongly suggests that the interactions between 
the individuals in the population occur through an active nutrient pool. For instance, in 
the case of interacting traders, a central exchange or reserve bank could represent such an 
active medium; this could selectively deplete or enrich individuals depending on taxation 
modes in the system. Our present study is based on such an example |5[; the detailed model 



of competing clusters on which it is based jo], forms the basis of our study here. The term 
'cluster' is used throughout this paper to refer a single coarse grained macroscopic unit that 
is a collection of numerous small microscopic subunits. The magnitude of such a cluster is 
characterized by its absolute mass. 

Beginning with the appearance of two seminal papers [3, Q], the study of complex net- 
works has attracted phenomenal attention from diverse fields. The spur in the growth of 
these studies began primarily with the prospect that numerous systems in nature can be 
compartmentalized into highly interconnected entities of dynamical units. To quote a few 
of the numerous examples that have a complex network topology : the Internet, social in- 
fractions, biological metabolic networks, predator-prey interactions in food chains, etc (see 
9|). Many real world networks, in spite of their inherent differences, have been found to 
have complex network topology, for example : social 10 1 and biological networks 



Predator-prey interactions seen in food webs and ecological networks are of special im- 



portance to our present study. O 
have a small world structure [l^. 



the numerous studies done on food webs, some of them 
13l |. while some other studies suggest that the food webs 
do not show either a small world or a scale-free topology but may have a different network 
structure jl^. While the exact form of the structure of a predator-prey system is subjective, 



the nature of the interaction has been shown to be highly nonhnear and therefore can be 
modelled mathematically 15 1. 



In the light of the above studies, it is becoming increasingly clear that interactions between 
agents are neither of the mean field type, where all of the species interact with all others, nor 
necessarily confined to local neighbourhoods. Therefore complex networks are a far more 
appropriate basis for the modelling of most systems with their combination of globality and 
locality, as well as their freer choice of agent-agent contacts. One of the main aims of this 
paper is thus to study the competitive cluster model of on networks. 



The plan of the paper is as follows. First, we introduce the model (Sec. [TTl); next, we do 
a deeper investigation of the problem in two dimensions in Sec. IIIII Here, we demonstrate 
the inherent complexity of the problem in finite-dimensional systems, where, contrary to the 
mean-field case, the largest cluster is not necessarily a survivor. Next, we investigate the 
main features of the model on a network. Starting with a regular network (which here is a 
finite-dimensional lattice) with nearest-neighbour interactions, we add non-local links to the 
existing lattice, with an associated probability p, in the range < p < 1 3]- We measure 
the survivor ratios as a function of the wiring probability p, as p is increased to reflect the 
topologies of small world networks and fully random networks {p = 1). Our system turns 
out to exhibit great sensitivity to small changes in the initial distribution of clusters; we 
characterize this using the concept of damage spreading (Sec. IV|l . 



In Sec. rvTl we ask the following question: can the destiny of a selected cluster be changed 
by networking? By 'networking' we mean there is a two way interaction between two or more 
clusters. We evolve a strategy to do this by selectively networking a cluster of our choice with 
non-local clusters based on their initial we masses. One of our most interesting observations 
in this paper is that a cluster which would die in its original neighbourhood, is indeed able to 
change its fate, becoming a survivor via selective networking. Next, in Sec. I VI 1 1 we calculate 
the pair-wise probability for a cluster to survive against larger and smaller neighbours; the 
probability distributions obtained show an interesting universality, independent of the mass 
considered. We discuss the implications of these and other results in our concluding section. 

Sec. Ivmi 
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II. MODEL 



The present model was first used in the context of cosmology to describe the mass ac- 
cretion of black holes in the presence of a radiation field [igI. Its applications, however, are 



considerably more general; used in the context of economics [5| , it manifested an interesting 
rich-get-richer behaviour. Here, we recapitulate some of its principal properties jsj. 

Consider an array of immobile clusters with time- dependent masses mi{t) located at 
equi-spaced regular lattice sites. The time evolution of masses of the clusters is given by the 
coupled deterministic first order equations. 

Here, the parameter a is called the mass accretion parameter and Qij defines the coupling 
between the clusters rrii and rrij; there is no self coupling. The first term in the R. H. S of 
Eqn. [1] represents the gain term, which is composed of two terms - the free rate of growth 
of an individual cluster (proportional to the parameter a > 1/2), and the rate of growth 
induced by all other clusters coupled through the surrounding medium. The coupling Qij 
between two clusters is proportional to the inverse square distance dij'^ between them at 
an initial time to 6|. The negative second term, — 1/mj, represents the dissipation to the 
surrounding medium. 

A logarithmic time is introduced in the study for convenience. We define a reduced time 
s = ln(t/to), where to is the initial time. Similarly, for convenience, we define reduced masses 
Xi = mi/t^^"^. Using the new variables, Eqn . [T]can be rewritten as. 



ds ' \ 2 ^"'^ V 2 

The primes denote the differentiation performed with respect to the reduced time variable 
s throughout. 

Continuing our recapitulation of the results of the model 6|, we consider a scenario 
where there is a single isolated cluster, whose initial reduced square mass is Xq. Under the 
dynamics defined by Eqn. [21 the cluster will be a survivor only if Xq > X^(= -y/^^j)' 
else it will eventually die out (see Fig. joj. Next, consider a system of two clusters with 
equal initial masses; it has been shown p] that there exists a critical coupling gc such that 
for g < they both survive if their initial masses are greater than X^. In the case of two 
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FIG. 1. (color online) The plots show the evolution of individual non- interacting clusters with a 
range of initial masses obeying Eqn. [2l Clusters with initial mass Xq greater than live forever, 



or else they die in time. Here a = 1.0, therefore X^, = ^/ ^ = -v/2 



clusters with unequal masses, the smallest of them dies first, in reduced time Si, and the 
largest cluster will either survive (if X{si) > X^) or die (if X{si) < X^). 

Consider the limit of an infinitely large number of clusters all connected to each other; 
this represents a limiting mean field regime, with fully collective behaviour involving long- 
range interactions. For g > gc, all the clusters but the largest will eventually die out. In 
the weak coupling regime {g < gc), which is the focus of this paper, the dynamics consist of 
two successive stages [6|. In Stage I, the clusters behave as if they were isolated; they grow 
(or die) quickly if their masses are greater (or less) than the threshold X^. In Stage II, slow 
collective dynamics leads to a scenario where, again, only the single largest cluster survives. 
This weakly interacting mean field regime shows the presence of two well-separated time 
scales, a characteristic feature of glassy systems Q]. 

Similar glassy dynamics also arise when the model is solved on a periodic lattice with 
nearest-neighbour interactions. The dynamical equations in Eqn. [2] take the form Gj: 




x„ 

by keeping the terms upto first order in g. Here m runs over the z nearest neighbours of the 
site n, where for a one- dimensional ring topology z = d, for a two-dimensional lattice z = 2d, 
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FIG. 2. The survival ratio S{s) plotted as a function of reduced time s, for clusters distributed in 
a regular one-dimensional lattice of size 100, 000. In Stage I, clusters grow independently, while in 
Stage II the growth is collective. Here, Si = 0.8 is the survival ratio at the end of Stage I and Soo 
is the asymptotic survival ratio. 



and so on for any space of dimension d. We summarize earlier results on the dynamics: In 
Stage I, clusters evolve independently; only those whose initial mass Xi{s = 0) > survive 
this stage. In Stage II, the dynamics are slow and collective; the survivors of this stage are 
fully isolated from each other by 'dead' sites, and their number asymptotes to a constant 
(Fig. [2]). These isolated clusters survive forever. 

Following the mean field scenario, where the largest cluster is the only one to survive, we 
ask if this would also be the case for the finite-dimensional case: that is, are the survivors 
the largest in the locality? Somewhat surprisingly, this is not always the case. While one 
can certainly rule out the survival of a cluster whose initial mass is less than threshold (X^), 
many-body interactions give rise to extremely complex dynamics in Stage II for clusters 
with X > X^. The asymptotic configurations of survivors in the finite-dimensional case are 
thus often nontrivial and counter-intuitive (see Sec. IIIip . 
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III. NONTRIVIAL DYNAMICS IN TWO-DIMENSIONAL REGULAR LATTICES 



Is 



In this section we do a deeper investigation of the cluster growth model [6|], with individual 
clusters placed on the vertices of a regular two-dimensional lattice with nearest-neighbour 
interactions (coordination number for all the clusters in the lattice is 2; = 4). While the 
results of Sec. |V]will show the extreme sensitivity of the model to the distribution of cluster 
masses, the results of this section will help us get a flavour of the complexity displayed by 
the model. 

We choose a 50 x 50 regular lattice for our investigations in this paper, unless otherwise 
specified; we have checked that there are no observable finite size effects for this choice. Initial 
masses are characterized by an exponential distribution P{Xq) = /iexp(— yuXo), where the 
mean is taken to be /i = — log(S'(i))/Xj,, for different values of ^i. The coupling parameter 
is set at g = 0.001 and the parameter a set to 1.0 jo] . 

First we consider two of the more obvious cases for demonstration (Fig. [3]). In Fig. [3](a), 
the central cluster has the largest mass among its neighbours, and in Fig. EJ^b), the central 
cluster has the smallest mass of those of its three neighbours which are larger than threshold 
(the fourth one will obviously not survive beyond Stage I). 

These two cases show results similar to the mean field scenario [gI; the central cluster 
with the largest initial mass becomes a survivor (Fig. [3]^a)) and the central cluster with the 
smallest initial mass dies off (Fig. El^b)). 

However, in general, such simplistic reasoning may not work. In a general scenario such 
as in Fig. |U consider the cluster located at co-ordinates (0,0), which has three larger and 
one smaller neighbours. Reasoning naively, we might expect that the larger neighbours 
X(i_o), 0), X(o,i) will eventually be able to kill the smaller central cluster X(o,o)- Con- 
trary to this, we find that in the asymptotic limit, the central cluster (0, 0) is a survivor. 
Similar results are obtained for the clusters at the coordinates (1, 1) and (1, —1); although 
surrounded by larger neighbours, they nevertheless survive against the odds. Another appar- 
ently counter-intuitive case is illustrated in Fig. |4]for the cluster with coordinates (—1,0). 
To start with, all the neighbours are smaller than the central cluster, but eventually the two 
neighbours at (0,0) and (—1, 1) survive at the expense of the central cluster, which dies. 

The explanation for such apparent anomalies is the many-body nature of this problem; 
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FIG. 3. A sample configuration of a subset of the system is shown here. The left panels show the 
cluster configurations at time s = 0. The right panels show the asymptotic configurations, (a) The 
central cluster dominates its immediate neighbours and emerges as the survivor, (b) Three larger 
neighbours eventually dominate the central cluster to become survivors by killing it. Open and 
filled circles represent non-survivors and survivors, respectively. 



the evolution of neighbouring sites is influenced by their neighbours, and so on. Examining 
the dynamical equation given in Eq. [3l we see that the overall contribution of the m = 4 
immediate neighbours of a given cluster at any instant of time is expressed by the summation, 

g dXmj , where Xm denotes the reduced masses of the m neighbours. This term 

should be positive for X'^ to be positive, which means that the overall rate of growth of its 
neighbours J2 ^'m should be negative. 

This can also be understood by inspecting Eqn. [2l In Eqn. [21 if ^ X'j is positive, then 
the overall growth rate of the neighbours of a cluster i under consideration is positive. Then 
the rate of growth of that cluster i, X'i, will eventually go negative. On the other hand, 
a negative ^ X'j will increase the rate X'i and enhance its growth. The precise manner 
in which the Yl -^'j evolves is complex given the many-body nature of the interactions (see 
Figs. [13] and [15] for sample time evolution plots of the combined rates of growth of the 
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FIG. 4. A small region of the two-dimensional lattice showing neighbouring clusters with com- 
parable initial masses. The numbers show the initial masses; open and filled circles represent 
non-survivors and survivors, respectively, in the asymptotic limit. This is an illustration of the 
complex dynamics referred to in the text. 

neighbouring clusters contrasted with the center cluster). Therefore, for a cluster to become 
a survivor, it should grow constantly; its rate of growth should thus always be positive. 

In other words: Clusters that survive stage I will be the ones with X > X^ and in stage I 
they grow independently; there is thus no reason for their rate to be negative. On the other 
hand in stage II, when there are interactions between clusters, even if a cluster temporarily 
has a negative growth rate, it will be fed on by its competitors and killed off. 

Next, we present plots showing the evolution of a given cluster and its immediate neigh- 
bours for the cases discussed in Fig. HI Consider the evolution of the cluster at (0,0), 
shown in Fig. |5] (a). The immediate neighbours of this cluster are seen to die when the 
central cluster grows (see Fig. EJ^b)). In Fig. M^c) we plot the rate of growth of the central 
cluster (in 'red') and the overall growth of the immediate neighbours SX'(jj) (in 'green'). 
Initially, both X'(o,o) and ,,•) are positive (Stage I), with the neighbours growing faster 

than the central site, in accord with intuitive expectations based on their initial masses. 
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FIG. 5. (Color online) (a) The central cluster -^(0,0) ^ ^ function of reduced time s. (b) The 
neighbouring clusters around (0, 0) as a function of s. (c) The rate of evolution of the central cluster 
and the combined rates of evolution of the immediate neighbours are compared - the central cluster 
wins the competition. 



When the collective dynamics of Stage II sets in, SX'(j j) falls and becomes negative; this 
is a demonstration of many-body interactions in action, where neighbours are influenced by 
their neighbours, and so on. Finally, soon after the neighbours begin to decay gets 
negative), the rate of growth of the central site, X'(o,o), shoots up, making it a survivor that 
continues to grow forever. 

Graphs similar to those in Fig. \5\ have been obtained for the other survivors of Fig. 
m viz. X(i 1) and X(i we have not plotted them here, in the interests of brevity. To 
complete our analysis, we show however the graphs in a contrasting case, where a given 
cluster eventually dies, despite being surrounded by smaller neighbours. Choosing X(i „i) in 
Fig. m we show the growth of the central site and its associated neighbours, respectively, in 
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FIG. 6. (Color online) (a) The growth of the central cluster g) as a function of reduced time 
s. (b) The growth of the neighbouring clusters around (—1,0) as a function of s. (c) The rate of 
evolution of the central cluster and the combined rates of evolution of the immediate neighbours 
are compared - the central site loses. 

Fig. El (a) and (b). Here, contrary to expectations, the central cluster dies in a finite time, 
while two of its neighbours manage to survive. These dynamics are seen clearly in Fig. [6] 
(c); in the fast dynamics of Stage I, the combined rate of growth of the neighbours is faster 
than the central cluster. After the onset of slow dynamics in Stage II, the rate of growth of 
the central cluster goes negative, which contributes to the growth of its neighbours. 

Clearly from the above, the simple rule of the mean field scenario, that the largest cluster 
is always a survivor, does not hold for lattice-based models with nearest-neighbour interac- 
tions. In fact, it is rare to be able to find cases like Figs. [3] (a)-(b), where we can ignore 
the effect of further neighbours. The generic case is much more as depicted in Fig. HI where 
information on the immediate neighbourhood of a site does not provide a reliable indication 
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of its fate, highlighting the innate complexity of our model on a lattice. 



IV. CLUSTER DYNAMICS IN COMPLEX NETWORKS 



In the previous section we discussed the complex dynamics exhibited by the cluster growth 
model on a lattice. However, both lattices and mean field scenarios are topologically unreal- 



istic - recent work indicates [121, Il3l . Il5| that most real- world systems are more appropriately 
solved on networks, which are spatially disordered arrays with non-trivial connections be- 
tween their sites. A particularly interesting example is the class of small world networks, 
which are constructed by starting with a regular lattice, and adding links randomly to its 
sites with probability p [7] . Alternative constructions are also possible, where existing links 
are 'rewired'; i.e. existing links are severed and then reconnected to randomly chosen lattice 
sites with probability p. Small world networks have the property that long-range interactions 
can co-exist with short-range interactions; they contain realistic elements like 'hubs', where 
certain sites are preferentially endowed with many connections, as is the case for example, 
in the aviation industry 

We shall in the following apply the former construction, whereby new links are added 
randomly with a given probability. In addition, we shall freeze the newly added links for 
all times, rather than evolving them continuously - which we call 'static wiring' 
result of this change, the average degree of the sites will be higher for all p > 0. 



18|. As a 



A. One-dimensional ring and two-dimensional square lattices 

First, we consider a regular one-dimensional ring lattice of size = 2000. To start with, 
the clusters located on the lattice sites evolve according to Eqn. [31 where the interactions 
are with nearest neighbours only. Next, we modify the lattice by adding new links between 
sites chosen randomly with an associated probability p G [0,1]. For p = 0, the network 
is undisturbed and therefore ordered, while for p = 1, the network becomes completely 
random. 

In our first scheme, we add links probabilistically starting with site i = 1 and end with 
i = N, only once: we call this the 1-cycle scheme. The survival ratios of clusters as a 
function of reduced time s for different values of wiring probability p are presented in Fig. 
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FIG. 7. (Color online) The plots of S{s) as a funwction of reduced time s, for increasing values 
of wiring probabilities p. The 1-cycle scheme is shown in (a) for a 1-dimensional ring and (c) for a 
2-dimensional square lattice. The 5-cycle scheme is shown in (b) for a 1-dimensional ring and (d) 
for a 2-dimensional square lattice. The insets in all the figures show the asymptotic survival ratios 
5*00 as a function of the probability p. Here, the system size for the 1-d ring is 2000 and for the 
2-d square lattice it is 50 x 50 - all our data is averaged over 10 random network configurations. 

[Tl^a). Consider the p = case, which corresponds to a regular lattice; here the survival ratio 
S{s) shows two stages. Stage I and Stage II, in its evolution. For all values of < p < 1, the 
existence of these well-separated Stages I and II is also observed. There is a noticeable fall 
in the survivor ratio as p is increased, though; this is clearly visible in the asymptotic values 
Soo{p) plotted with respect to p in the inset of Fig. [7](a). As the probability p increases, the 
number of links increases leading to more interaction and competition between the clusters 
and hence a decrease in the number of survivors. As expected, this behaviour interpolates 
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FIG. 8. The plot of asymptotic survival ratios Soq{s) as a function of the probabilities p for 1- and 
5-cycle schemes in 1 dimension and 2 dimensions. The error bars for all the graphs do not exceed 
0.003 and are smaller than the plot symbols. 

between the two characteristic behaviours relevant to the regular lattice and mean field 
scenarios. 

Next we implement the 5-cycle scheme, where the rewiring is done five times. Figure 
Wljo) shows the survivor ratio S{s) as a function of s for this scheme. We can clearly see the 
decrease in the number of survivors at all times, for increasing values of p. The asymptotic 
survival ratios 6*00 for the 1-cycle scheme and 5-cycle scheme are shown in the insets of Fig. 
Wi^) and (b), respectively, where the monotonic decrease of S^o with increasing p is very 
clear for both cases. In addition, for all values of p, the survivor ratio in the 5-cycle scheme 
is consistently smaller than in the 1-cycle scheme. 

The same procedure for a two-dimensional square lattice of size 50 x 50 is followed, and 
survival ratios obtained as a function of s using the 1-cycle scheme and 5-cycle scheme 
(see Figs. [71^c)-(d)). The asymptotic survivor ratios for these two cases are shown in the 
insets of Figs. [71I^c)-(d); they follow a decreasing trend with increasing p, similar to the one- 
dimensional case. We plot the asymptotic survival ratios of the 1-cycle and 5-cycle schemes 
in 1-dimensional and 2-dimensional lattices in Fig. [8] for comparison. Notice that the trends 
in the plots for the 5-cycle scheme are far smoother than those for the corresponding 1-cycle 
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plots. This is expected because of better statistics achieved in the 5-cycle scheme. 



V. DAMAGE SPREADING RESULTS 



In this section, we measure the sensitivity of the system to shght perturbations in terms 



of damage spreading, 
of biological systems 
mechanical systems 
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'he concept of damage spreading was introduced in the context 
and was widely used in the study of the dynamics of statistical 



2l| . A system is said to have damage spreading if the Hamming 



distance between two replicas, dictated by the same evolution equations, under slightly 
different initial conditions, grow with time. This concept has also been used to study phase 
transitions in random networks [22] . However, our definition of damage spreading is different 
from the usual definition based on the Hamming distance, used by workers on Ising models; 
we will elaborate on this below. 



A. Sensitivity to small perturbations 

Again, we consider a 1-dimensional ring lattice with only nearest-neighbour interactions. 
An exponential distribution of initial masses of clusters is put on the lattice, as before. 
Starting with two identical initial mass configurations Qa and fi^, we swap two randomly 
chosen clusters in the configuration Qb- We recall from the discussion above that the 
evolution of clusters in a finite-dimensional lattice consists of two stages, fast and slow, and 
that an initial cluster can be a survivor at the end of Stage I only if its mass Xi > X^. 
Therefore, while perturbing the initial mass configuration Qb, we consider a variety of 
regimes that takes into account all possible combinations of clusters that are swapped. 
Table [T] summarizes all such combinations of X^ and Xj, scaled with respect to the threshold 
square mass X^. Next, we evolve the two configurations Qa and fi^ in time and find the 
absolute difference of the total mass of the two configurations, |m^ — m^l, with time. Figure 
|9] shows the difference \mA — TnB\ for all five regimes shown in table HI Note that irrespective 
of the regime, the quantity — m^l diverges with time. We conclude therefore that the 
system is very sensitive to the initial masses of the clusters, even for the regular lattice, 
where no random links are added. 

We now summarize the above results. First, even for the regular square lattice, we have 
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Regime Xi 



Regime 1 « X^, 
Regime 2 < X^, 
Regime 3 < 
Regime 4 > 



«X, 
«X^ 

>x^ 
>x^ 



Regime 5 > > X^,{Xi ~ Xj 



TABLE I. The five regimes defined liere are with respect to the threshold mass X^,. The clusters 
located at sites i and j are swapped in the configuration Q b ■ 

demonstrated the complexity of the dynamics: it is impossible to predict the fate of a 
randomly chosen cluster unless its initial mass is less than the critical mass X^. In Sec. IIVI 
we have started our study of the model on a network structure by adding new links (static 
wiring) the square lattice; we have shown that as the probability of random connections 
p increases, the asymptotic number of survivors decreases monotonically. This decrease, 
attributable to the death of clusters with marginal masses when their contacts increase in 
number, matches with the results of the mean field scenario js]. From a damage spreading 




FIG. 9. (Color online) The absolute mass difference \mA — msl of two mass configurations Qa and 
flB, as a function of the reduced time s. The case shown here is for a regular 1-dimensional regular 
ring lattice of size 1000. 
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point of view, we also see that the system displays great sensitivity to the distribution of 
mass within it; the average mass difference between a perturbed and parent cluster diverges 
in time, irrespective of the masses of clusters that are swapped. 

Given all of this, we might guess that it would be even more difficult to predict the future 
of a randomly chosen cluster (again, for Xj > X^) when the model is solved on a network. 
In the next section, we attempt this in order to address the following question - can selective 
networking be used to change the fate of a given cluster? 

VI. NETWORKING STRATEGIES 

Since the central issue in this model is the survival of clusters against the competition, 
the most interesting use of networks would be in the finding of a smart networking strategy 
which is able to change the fate of a cluster. This can be accomplished by making a marginal 
survivor into a strong one, and more dramatically, by seeing if a dying cluster can be revived 
to life. 

We systematically investigate the effect of adding a finite number of nonlocal connections 
to a chosen central cluster. In Sec. IIIH it was shown that the growth or decay of a cluster 
is solely dictated by its relative rate of change versus the cumulative rate of change of 
its neighbours. The key to better survival should therefore lie in choosing appropriately 
decaying non-local clusters to network with; we will show in the following, that the fate of 
a cluster can indeed be modified if this is done. 

Given the complexity of the problem, we cannot categorize the clusters based on their 
rates of growth; the most viable option is a categorization of the clusters based on their 
initial masses. We divide the nonlocal connections into two classes: viz., class A comprises 
eventual non-survivors {X < X^), while class B comprises would-be survivors {X > A^) ). 
In the next subsection, we look at the outcome of networking with members of class A. 

A. Nonlocal connections with eventual non-survivors (Aj < A^,) 

Recall that non-survivors (< A*) die very early during Stage I. In connecting such clusters 
to a given cluster with A > A^, we can be sure that they will never be able to compete with 
it, much less kill it. Could such 'harmless' clusters, nevertheless be useful in enhancing the 
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FIG. 10. (Color online) (a) The central cluster here is a survivor in its original configuration, 
(b) The asymptotic mass of the central cluster diverges with increasing connectivity with eventual 
non-survivors. 

survival chances of clusters with X > X^l 

Let us consider a central cluster that will be an eventual survivor on its own in a given 
neighbourhood; one such case is shown in Fig. [TDl We now let the central cluster network 
with eventual non-survivors from all over the lattice, and record its growth as a function of 
the number of the clusters it networks with; the results are shown in Fig. [10] (b). When the 
neighbours of the central cluster die, their contribution goes to zero (see Eqn. 3) and the 
solution to the resulting first order linear differential is therefore an exponential as shown in 
Fig. [10] (b). It is evident that there is a marked increase in the mass of the central cluster 
as more small clusters are connected to it; since none of these will survive past Stage I, 
their cumulative growth rate will be negative, leading to an increase in the rate of growth of 
the chosen cluster. This makes it a better (more massive) survivor asymptotically. As the 
number of nonsurviving contacts is increased, the mass of the chosen cluster increases (Fig. 

We now use this observation to try to change the fate of a dying cluster, and turn it 
into a survivor. Although the nonlinear, many-body nature of the problem makes it rather 
difficult to construct such cases, we have been able to network an erstwhile dying cluster 
with suitable nonlocal contacts, and returning it to life. Figure [TT] shows our results: the 
central cluster would eventually have died in its original environment, but on adding 100 
small clusters (whose mass < X^) it manages to live. Further additions, e.g. 200 or 364 
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clusters evidently make it a better survivor (Fig. [TT] (b)). 



B. Networking with would-be survivors Xi > 

Here we are dealing with a far more complex problem than the case discussed in earlier 
subsection. When we choose clusters with masses larger than X^,, their lifespan will exceed 
Stage I; thus during Stage I they will have a positive growth rate. Depending on their 
individual environments, they may survive through Stage II - thereby having a positive 
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FIG. 11. (Color online) (a) The central cluster dies in its original configuration without extra 
connections. On addition of more and more clusters it becomes a survivor, (b) A crossover is seen 
here as the central cluster becomes a survivor from being a nonsurvivor, with the increased number 
of connections to non-surviving clusters (< X^). 
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FIG. 12. The central cluster is networked with would-be clusters (Xi > X^) that are less massive 
than the central cluster (Xcm)- (a) The initial configuration, where the central cluster has 6 
nonlocal connections, (b) the asymptotic state with the only survivor being the central cluster, (c) 
With the addition of one more nonlocal cluster to the existing configuration (d) the central cluster 
dies. In (b) and (d), the open circle represent nonsurvivors and dark circles represent asymptotic 
survivors. 

growth - or die as a result of a negative growth rate. Networking with such clusters is 
therefore a rather delicate matter. 

First, we consider nonlocal connections with would-be survivors (Xj > X^) which are less 
massive than our chosen cluster (Xj < Xcm). Such would-be survivors will live beyond Stage 
I and will grow with a positive rate. This is shown clearly in Fig. [T] where clusters with 
Xj < X^ die in finite time and clusters with Xj > grow beyond stage I. From a mean 
field perspective, we would therefore expect to see a decrease in the chances of survival of 
the chosen cluster as it networks with more and more would-be survivors. 

Figure [12] shows a sample scenario, where the central cluster is connected with nonlocal 
clusters which are not only would-be survivors (Xj > X^,) , but also less massive than the 
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central cluster. In Fig. [12] (a) the central cluster networks with 6 would-be survivors and 
is able to survive in the asymptotic state (Fig. [12] (b) ). On the other hand, adding one 
more would-be survivor (Fig. [12] (c)) to the existing network of the central cluster kills it in 
the asymptotic state (Fig. [T2] (d)). We also note that, this change in the fate of the central 
cluster - from being a survivor to being a dying cluster, due to the additional connection - 
also changes the fate some of the clusters linked to the central cluster (contrast Figs. [12] (b) 
and (d)). 

To understand the dynamics leading to the observed behaviour we look at the rates of 
growth of the central cluster and its neighbours. A sample result corresponding to the cases 
shown in Fig. [T2]is discussed in Fig. [12] Here, we see that increasing the number of nonlocal 
connections with would-be survivors leads to a fall in the absolute value of X™ as well as 
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FIG. 13. (Color online)The central cluster of a configuration given in Fig. [121 a) is networked 
with clusters whose Xn > X-^. In (a) the growth of Xcm with an increasing number of connections 
is plotted - n = 4 corresponds to a regular lattice. There is a crossover seen when the number 
of connections increases from n = 10 to n = 11; for larger n values, the central site dies. This 
observation is supported by the rates of growth of Xcm and its neighbours X„ (b) when n = 10 
and (c) when n = 11. 
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2 nonlocal connections 
(a) Initial configuration (b) Final configuration 




3 nonlocal connections 




(c) Initial configuration (d) Final configuration 

FIG. 14. The central cluster is networked with would-be clusters {Xi > X^) that are more massive 
than the central cluster [Xcm)- (a) The initial configuration with 2 nonlocal connections, (b) and 
its asymptotic state; the only survivor here is the central cluster, (c) With the addition of one 
more nonlocal cluster to the existing configuration (d) the central cluster dies. In (b) and (d), the 
open circle represent nonsurvivors and dark circles represent asymptotic survivors. 

its rate of growth X'^^. Beyond a certain number of networked contacts, the central cluster 
begins to decay, and eventually dies. This crossover from life to death happens when the 
cumulative rate of growth of the neighbours is larger than that of the central cluster 

X^^. Unfortunately the intricate many-body nature of this problem precludes a prediction 
of when such crossovers might occur in general. We will have more to say on this issue in 
the next section. 

Finally, in the case where a given cluster networks with larger would-be survivors {Xi > 
X^ and Xi > Xcm), one would expect a speedier death. One such sample scenario is depicted 
in Fig. [THand the corresponding rates of evolution of the clusters in Fig. [TBI We notice that 
in its original configuration the central cluster (with four neighbours, (n = 4)) is a survivor. 
As we increase the number of networked connections, its growth gets stunted; there is a 
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substantial fall for 2 extra links (n = 6 in Fig. [T5|) . Adding one more link (n = 7) does the 
final damage; the central cluster dies. The rates shown in Fig. [15] (b) and (c) for = 6 
connections and n = 7 connections vividly capture the competition for survival, leading to 
life in one case and death in the other. 

As expected, we observe that fewer connections (here, n = 7) are needed, compared 
to the earlier case with smaller would-be survivors {n = 11), to kill the chosen cluster. 
In closing, we should of course emphasise that these results are sample results chosen to 
illustrate a qualitative point: the choice of different nonlocal sites to network with would 
certainly modify our n values. 
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FIG. 15. (Color online) The central cluster of a configuration given in Fig. [T^b) is connected 
with nonlocal clusters with > Xcm- lu (a) the growth of Xcm with an increasing number of 
connections is plotted - n = 4 corresponds to a regular lattice. There is a crossover seen when the 
number of connections increases from n = 6 to n = 7; for larger n values, the central site dies. This 
observation is supported by the rates of growth of Xcm and its neighbours X„ (b) when n = 6 and 
(c) when n = 7. 
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VII. UNIVERSALITY IN THE DISTRIBUTION OF SURVIVORS - A CLUE TO 
SURVIVAL AGAINST THE ODDS? 



The results of the previous section indicate that if a cluster wants to change its fate, 
it should always network with non-survivors; however, it seems intuitively reasonable to 
suppose that it might also be able to do this if it networks with would-be survivors, provided 
that they are much smaller than itself. Finding such cases is possible ex post facto; their 
prediction is next to impossible given the complexity of the problem. 

An interesting case to consider might be one where, against the odds, the most massive 
cluster in a given neighbourhood dies marginally, because of the unexpected growth of its 
erstwhile smaller neighbours. If, now, it is able to network with potential survivors much 
smaller than itself (who, as a result of their own dynamics, would eventually perish in their 
own neighbourhoods), this would be a promising way of selectively networking itself to life. 
Similar issues are associated with survivors who survive against the odds; they could be 
killed off by suitable networking. 

A systematic way of exploring this question is to look at the statistics associated with 
some of these rare events. We look first at the four immediate neighbours of a given cluster, 
and consider their pairwise interactions with it. Clearly, had such a pair been isolated, the 



larger cluster would have won jo]. However, many-body interactions in the lattice mean that 
this is not always true. We therefore ask the question: what is the proportion of cases where 
the larger cluster wins? 

Each survivor has four neighbours; we first calculate the probability distribution of the 
initial mass differences in a pairwise fashion between a survivor and each of its neighbours. 
The initial mass differences are given by 6Xi = Xcm — Xi {i = 1, 2, 3, 4) corresponding to the 
four neighbours - right, left, bottom and top - of a survivor. The distribution of SXi for all 
the survivors is shown in Fig. [161 Here, a negative 6Xi means that the survivor is smaller 
than its neighbour, and conversely for positive 6Xi. All four distributions corresponding 
to four neighbouring pairs overlap due to isotropy; the resulting distributions are universal 
functions of mass differences, depending only on . 

We also obtain the cumulative mass difference between survivors and all of their four 

neighbours viz. AXcm — X]^=i ^« = Ylt=i ^-^i (^^^ -^i§- [13) • The distributions of Ylt=i ^« 
plotted in Fig. [17] for different values of ii. For a positive cumulative mass difference we 
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FIG. 16. The plots show the distribution of pairwise mass differences between survivors and their 
four neighbours. They are obtained for exponential distributions of initial masses with different 
mean values where /i = — log(5'i)/X^. The plots are for (a) = 3.92 (5i = 0.6) (b) 5.607 
{Si = 0.7) (c) 8.963 {Si = 0.8) and (d) 18.982 (5i = 0.9). The system size is 400 x 400. 

know that the survivor has a larger initial mass than the sum of its neighbours, matching 
our intuition based on the mean-field regime. The negative side of the distribution is more 
interesting, comprising survivors whose initial masses are smaller than the sum of the initial 
masses of their immediate neighbours. 

Notice that both the survivor-neighbour pair distribution (Fig. [T6l) . and the survivor- all 
neighbours cumulative distribution (Fig. [17]) get broader with increasing /z. This is because 



increasing /i = — log(5'i)/X^ 6<] increases the number of potential survivors 5*1 beyond Stage 
1. In each case, the fraction of area under the negative side of the survivor pair-distribution 
gives an estimate of survivors against the odds - an example of some of the rare events 
alluded to at the beginning of this section. 

Figure [18] shows this fraction, both in terms of individual survivor-pair distribution and 
cumulative distribution, as a function of the /i of the initial mass distribution, for different 
system sizes. There are more survivors, hence more survivors against the odds, leading to 
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FIG. 17. (Color online) The plots show the distribution of mass differences between survivors 
and all of their four neighbours. They are obtained for exponential distributions of initial clusters 
with different mean values l//i, where fi = — log(S'i)/X^. The plot in 'red' represents the case for 
1//X = 3.92 {Si = 0.6), 'green' represents l/fj, = 5.607 {Si = 0.7), 'blue' represents l//i = 8.963 
{Si = 0.8), and 'pink' represents l//i = 18.982 {Si = 0.9). The system size is 400 x 400. 



an increase in the fraction plotted on the y-axis of Fig. [18] for both distributions. For the 
largest system size, there is full isotropy in the pairwise distributions; the probability of 
finding a survivor against the odds is now seen to be a regular and universal function of 
H in both pairwise and cumulative cases, relying only on mass differences rather than on 
masses. Finally, the cumulative distribution gives a more stringent survival criterion than 
the pairwise one, as is to be expected from the global nature of the dynamics. 



A major conclusion to be drawn from Fig. [18] is the following: there are clusters which die 
against the odds (clusters which are more massive than the survivor in either a pairwise or 
cumulative sense). Recalling again that the distributions plotted are independent of mass, 
it should therefore be possible to pick dying clusters whose mass was greater than threshold 
and to network them with their peers so that they can survive. This will be the subject of 
future work. 
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FIG. 18. The plots show the probability of finding a successful central cluster with a larger 
neighbouring cluster, to start with. The probability increases with increase in Si (refer to Figs. [T6l 
and I17p . While the solid line show the combined cumulative probabilities (CCP) (from Fig. [T7|) . 
the rest of the lines represent the individual cumulative probabilities (from Fig. [T6|) . 

VIII. DISCUSSION AND CONCLUSION 

In this paper we have investigated the cluster growth model of [sl on a square lattice. 
Our main emphasis has been on stochastic networking strategies: long-range connections 
are introduced with probability < p < 1 to existing lattice sites. 

We find first of all that the qualitative features of the networked system remain the same 
as that of the regular case, in that the presence of two well-separated dynamical stages is 
retained; the glassy dynamics and metastable states of the earlier work jo] persist. However, 
the number of survivors decreases as expected with increasing p, in accordance with a mean 
field perspective. Also, the system is very sensitive to the smallest changes in the initial mass 
configuration; the divergence is, interestingly, exponential irrespective of the magnitudes of 
the masses swapped. To quantify this, we have used the concept of damage spreading. 

The central result of this paper is the use of smart networking strategies to modify the 
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fate of an arbitrary cluster. We find that it is safest to network with non-survivors; their 
decay and eventual death lead to the transformation of the destiny of a given site, from 
death to survival, or from survival to stronger survival. Networking with peers or with those 
of larger initial mass in general leads to weakening, and an almost inevitable death. 

However, the above is not immutable: the probability distributions in the last section of 
the paper indicate an interesting universality of survival 'against the odds'. In future work, 
we will exploit this result to develop schemes for survival strategies; in complex systems of 
this type, it seems clear that the phrase 'survival of the fittest' is more hkely to be replaced 
by 'survival of the best networked'. 
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